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We show that the pseudospin being an additional degree of freedom for carriers in graphene can 
be efficiently controlled by means of the electron-electron interactions which, in turn, can be manip- 
ulated by changing the substrate. In particular, an out-of-plane pseudospin component can occur 
leading to a zero-field Hall current as well as to polarization-sensitive interband optical absorption. 
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Introduction. The charge carriers in graphene are de- 
scribed at low energies by an effective Hamiltonian be- 
ing formally equivalent to the massless two-dimensional 
Dirac Hamiltonian Hq = hvo{h'axkx +cryky), where 
!/ = ± refers to the two inequivalent corners K, K' of 
the first Brillouin zone, vq « lO^ms^^ is the effective 
"speed of light", k is the two-component particle mo- 
mentum operator, and a^^y are Pauli matrices describ- 
ing the sublattice degree of freedom also referred to as 
the pseudospin In the original Dirac Hamiltonian the 
Pauli vector a represents the spin of a spin- 1/2 particle 
which can be detected in Stern-Gerlach-like experiments. 
The pseudospin in graphene is formally similar to the true 
electron spin with an important distinction given by the 
behavior under time and parity [3| inversion. For the 
above effective model the time reversal operator T is just 
the operator C of complex conjugation, T — C, fulfilling 
THqT^^ — Hq'^, and the operators Hq and Hq get in- 
terchanged. On the other hand, if one would (formally!) 
interpret the Pauli matrices as components of a genuine 
spin[5| (i.e. an angular momentum), the time reversal 
operator would read f = ayC giving TH^f-^ ^ H^. 
The parity operator P flips the sign of the spatial co- 
ordinates interchanging the two sublattices and, similar 
to T, fulfills PH^P-^ = H-". Thus, the initial Hamil- 



tonian Hq 



Hq is PT-invariant but, as we shall see. 



the exchange interactions can break either invariance. It 
also is important that the pseudospin is not linked with 
the internal magnetic moment of an electron and does 
not directly interact with the external magnetic field pro- 
hibiting Stern-Gerlach type experiments. In contrast to 
that, we predict situations where the pseudospin mani- 
fests itself in observable quantities and can be detected in 
transport as well as optical measurements on graphene. 

First of all we show that the exchange electron-electron 
interaction can alter the pseudospin orientation in a very 
broad range. In an eigenstate of Hq the pseudospin is 
always in the xy-pl&ne. As we shall see shortly, the ex- 
change interactions can turn the pseudospin texture to 
the out-of-plane phase with the out-of-plane angle de- 
pending on the absolute value of the particle momen- 
tum. This is due to the huge negative contribution to the 
Hartree-Fock ground state energy from the valence band 
(i. e. "antiparticle" states) which cannot be neglected 
in graphene because of the zero gap (i. e. zero effective 



mass of carriers) and large effective fine structure con- 
stant a* = e'^/{ehvo) where e is the dielectric constant 
depending on the environment Q . The exchange contri- 
bution to the ground state energy has previously been 
studied in both monolayer and bilayer graphene regard- 
ing properties such as the electronic compressibility 
and ferromagnetismjsl-fiol , but the importance of the in- 
terplay between pseudospin and electron-electron inter- 
actions has been recognized only in ll| where single layer 
graphene was mentioned in passing. 

Having established the possibility to create an out-of- 
plane pseudospin orientation by means of the exchange 
interaction, we apply the Boltzmann approach to derive 
the electrical conductivity tensor which turns out to have 
Hall components even though the external magnetic field 
is absent. The mechanism of this phenomenon is in- 
timately linked to the pseudospin-momentum coupling 
which can be read out immediately from the Hamiltonian 
Hq . Similar to the skew scattering of electrons on impu- 
rities in spin-orbit coupled systems partly responsible for 
the anomalous Hall effect 12 , l3| , the carriers in graphene 
do also skew to one side of the Hall bar as long as their 
pseudospin has non-zero out-of-plane component. This 
effect has been intensively studied ^l4fil6i] assuming that 
the out-of-plane component occurs due to the band gap 
opened by spin-orbit coupling IJ] which, however, seems 
to be weak in graphene [2|. We emphasize that neither 
spin-orbit coupling nor an external magnetic field is nec- 
essary to obtain a Hall current in graphene being in the 
pseudospin out-of-plane phase. 

Experimental manifestations of the pseudospin are not 
limited to the electron skew scattering phenomenon but 
can also be seen in the interband optical absorption. Per- 
forming optical measurements on graphene [l7| one can 
obtain direct information regarding conduction and va- 
lence band states without advanced sample processing 
necessary for transport investigations. Moreover, the pe- 
culiar properties discovered so far make graphene a very 
promising material for optoelectronic applications [l8|. 
Optical absorption via the direct interband optical tran- 
sitions in graphene has been investigated in [l^ but 
the mechanism considered there lies essentially in the 
two-dimensional nature and gapless electronic spectrum 
and does nor directly involve the pseudospin orientation. 
Here we show that, due to the out-of-plane pseudospin 
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orientation, the interband absorption can be substan- 
tially reduced or enhanced as compared to its universal 
value Tre^/fi-c just by switching the helicity of the circu- 
larly polarized light. 

Exchange interactions. The Coulomb exchange Hamil- 
tonian is given by 

with ?7[k-k'| — 27re-^/e|k — k'l and k' = ± being the band 
index with k = + for the conduction band. The inter- 
valley overlap is assumed to be negligible, and the eigen- 
states ofH" = H^+H^^^^ can be formulated as ^'^^(r) = 
e^'^lxife) with spinors |x+fe) = (cos ^, z/sin ^e''''^)^, 
Ix'Lk) = (sin ^,-1/ cos ^e''^'^)'^, and tan(p = ky/k^. 
Thus, a non-zero out-of-plane pseudospin component cor- 
responds to 7^ 7r/2. To diagonalize 77" the following 
j'-independent equation for iSk must be satisfied (20j| 

/cP k' 
-^--jK't/|k-k'| [cos'&k' sinz9fc- 

K' 

— siniJfc' cosi^fc cos{ip' — ip)] = 0, (2) 
where the integration goes over the occupied states. Note 



that the conduction and valence states are entangled, 
and the latter cannot be disregarded even at positive 
Fermi energies. Thus, in order to evaluate the integrals 
in Eq. ([2]) a momentum cut-off A is necessary. Its value 
~ O.lnm"^ is usually chosen to keep the number of states 
in the Brillouin zone fixedQ, but our outcomes do not de- 
pend on any particular choice of A. Substituting x = fc/A 
we arrive at 

Attx cos 7?t , , 

= (3) 

a* ^ ' 

2-ix 1 

/. [ , . COS sin i^^ ~ sin iJfe/ cos i}}^ cos kp' 

dip ax X . 
J W + x'^ — 2xx' cos ip' 

kp/A 

The momentum cuf-off is obviously much larger than the 
Fermi momentum fc^? at any reasonable electron doping, 
and therefore we can set the lower integral limit to zero. 
Besides a trivial solution with = tt/2 independent of 
k, there are non-trivial ones i?i = i9(fc) and i?2 — TT — -d{k) 
with i^(fc) shown in Fig. [T] for different a*. The solutions 
i9o and i9i.2 represent to two phases with different total 
ground state energies iSJ^t {E^^^) for the in-plane (out-of- 
plane) pseudospin phase. The difference AE^tot — -Etot ~ 
E°^^ per volume for a given spin and valley reads 



A^tot 
HvqA^ 



27r ^ ' 



dip / dtp' dx dx'a 



I (1 — siniJfe' sint?fe) cos{p' — p) — cos-dk' cosi?fc 



327r3y^ 



2xx' cos((^ — p') 



(4) 



The energy difference for a* ^ 1 is small because the 
integrand in Eq. Q is always multiplied by x' and there- 
fore vanishes at x' — 0, but at larger x' the -dk' gets close 
to 7r/2, and the integrand vanishes again. The Inset in 
Fig.[T]shows, however, that strong electron-electron inter- 
actions make the out-of-plane phase energetically prefer- 
able. The estimates of a* for clean graphene vary from 
2 (Ref. (il) to 2.8 (Ref. [3) and are on the borderline of 
the out-of-plane phase. Moreover, the presence of dis- 
order can change this qualitative picture essentially Q. 
Most importantly, Eq. ^ is valid for both valleys and 
both solutions 'i?i^2- Thus, it is possible to choose either 
the same or opposite solutions for two valleys. The for- 
mer choice breaks the parity invariance whereas the latter 
one does so with the time reversal symmetry. Both cases 
are worthy of consideration. 



The single-particle spectrum is independent of the val- 



ley index and given by 

2tt 1 

Ek{x) . ^* f ^1 f f ^ ' ' 

— — — = Kx sm Vk — ■:— I dp I dx x 



1 — K(cosi9fc/ cos?9/£ + sin-iJfc/ sini^fc coslys') 

X , = , (5) 

y + a;'^ — 2xx' cos p' 

and the group velocity can be written as = 
Wk(cos (y9, sin lys)-^ with being 



— = Ksm'dk + — / dp' / (ia::'a::'a;(l — cos (js') 
vo An J ^ J 



1 — K(cosi?fc' cos-dk + sinzJfc/ aindk cos 09') 

X ^ 3 -. (6) 

{x'^ + x'^ — 2xx' cos p') 2 

The dispersion law ^ is depicted in Fig. [2] for graphene 
placed on Si02 substrate. The interactions shift the 
bands down to lower energies and change the density 
of states but, most importantly, they open a gap 21 1 
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FIG. 1: The pseudospin out-of-plane angle ^{k) for differ- 
ent environments numerically calculated from Eq. @- The 
corresponding values of the substrate-dependent effective fine 
structure constant a* are taken from Ref. Q|. The inset shows 
the total ground state energy difference (|4| between the in- 
plane and out-of-plane phases for different effective fine struc- 
ture constant a* = /ehvo- Increasing a* makes the out-of- 



plane phase more preferable. 
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FIG. 2: The dispersion law Enik) in the in-plane (dashed 
curves) and out-of-plane (solid curves) phases for a* = 0.8 
corresponding to Si02 substrateQ. The curves for both 
phases coincide for momenta larger than a certain critical 
value where "dk = 7i"/2 becomes independent of fc, see Fig. [1] 
The inset shows the gap region in detail. 



between the valence and conduction band as soon as the 
system changes to the pseudospin out-of-plane one phase. 

2 1 

The gap at fc = equals ^-^ / dx' cos d]^' . Note that the 

^ 

group velocity ^ vanishes at small momentum fc/A <C 1 
as long as the system is in the out-of-plane phase corre- 
sponding to the almost flat bands close to fc = shown in 
the inset of Fig. [2] From now on we assume n-doping so 
that the Fermi energy is always higher than the bottom 
of the conduction band. 

Zero-field Hall current. To describe the Hall conduc- 
tivity due to skew scattering we utilize the semiclassical 
Boltzmann approach which allows a physically transpar- 



ent interpretation of this mechanism isl 15 1. In general 
the anomalous Hall conductivity contributions can be 
classified by their mechanism: (i) The intrinsic contri- 
bution is due to the anomalous velocity of carriers (be- 
ing non-diagonal with respect to the band index |22| ) 
which is coupled to the equilibrium part of the distri- 
bution function, (ii) The side-jump contribution follows 
from coordinate shifts during scattering events. It occurs 
in the non-equilibrium part of the distribution function 
as well as in the anomalous velocity [H [H. (iii) The 
skew scattering contribution is independent of the co- 
ordinate shift and of the anomalous velocity. It occurs 
when the scattering rate is asymmetric with respect to 
the initial and final states and, therefore, must be con- 
sidered beyond the first Born approximation The first 
two conductivity contributions do not depend on disor- 
der but on the out-of-plane angle -dk and can be adopted 
from [l^. Here, we focus on the skew scattering contri- 
bution which can be described using the interband inco- 
herent Boltzmann equation where the anomalous velocity 
is neglected but the scattering probability is calculated 
up to the third order in the short-range scattering po- 
tential with the momentum-independent Fourier trans- 
form V. In linear order in the homogeneous electric field 
E this equation reads -eEv^ [-df°{Ek)/dEk] = ^[/fc], 
where f'^{Ej.) is the Fermi-Dirac function, is the non- 
equilibrium addition, and v^, E^ are given by Eqs. (|6l5p 
with K = +. The collision integral can be written as 
^[/k] = / f|:;^w'kk' (/k' - /k) with Wkk' being the scat- 
tering probability. We divide Wkk' into two parts. The 
first one is proportional to the cosine of the scattering 
angle and calculated up to the second order in V. The 
second one is proportional to the sine of the scattering 
angle and calculated up to the third order in V. These 
two parts correspond to the conventional and skew scat- 
tering respectively which can be alternatively expressed 
in terms of the momentum relaxation times, cf. Ref. [3] 

(r|f)-i = n,fcy2(i + 3cos2,?^)/(4ft\fe), 
(tI)-^ = vmk^V^ cos dlsin^ dl/{?>h^vl). (7) 

Here, Ui is the concentration of such scatterers. Since 
oc whereas t^^ oc l/V"^ it is natural to assume 

^ Tip and the Hall conductivity for a given valley 



can be estimated as ~ o-xxt!^ /t'^ \k=kp which can 
vary in a quite broad but finite range because neither 
of r's diverges at low doping thanks to the fc-dependent 
group velocity Q . Note that K and K' contribute iden- 
tically to the total Hall conductivity ayx = '^yx if 
the out-of-plane pseudospin polarization is opposite in 
the two valleys, i. e. and are assigned to iJj!' and 
respectively, and the time reversal invariance is bro- 
ken by the exchange interactions. On the other hand, if 
the out-of-plane pseudospin polarization is the same in 
both valleys (i. e. either of i?i_2 is assigned to both 
breaking the parity invariance) the Hall currents in the 
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two valleys have opposite directions resulting in the val- 
ley Hall effect |23p — another analog of the well known 
spin Hall effect [14]. 

Interband optieal absorption. From Hq one can de- 
duce the following interaction Hamiltonian between the 
electromagnetic wave and carriers in graphene = 
— {vaxAx -\-(JyAy) which couples the vector potential A 
and pseudospin a. As consequence, the inter-band tran- 
sition matrix elements turn out to be sensitive to the light 
polarization and pseudospin orientations in the initial 
and final states. To be specific we assume monochromatic 
light of frequency uj, normal incidence (i.e. zero momen- 
tum transfer from photons to electrons) , and circular po- 
larization (fulfiUing Ax = ±iA/y/2, Ay = A/y/2). The 
probability to excite an electron from the valence band to 
an unoccupied state in the conduction band can be cal- 
culated using the golden-rule. Finally, the absorption P" 
can be calculated as a ratio between the total electromag- 
netic power Wa absorbed by graphene per unit square 
and the incident energy fiux Wi ~ uP A^ / Aire. Then, the 
optical absorption for K valley {u = +) reads 



P+ 



7re^ 4Awo 
he (jj 



dxx ■ 



sm 



4 K 



it 



hvoA 



(8) 

where the multipliers sin'*(i?^/2) and cos'*(i?^/2) are for 
two opposite helicities of light, and for K'-valley they are 
interchanged. If the out-of-plane pseudospin polarization 
is chosen to be opposite in the two valleys, then the to- 
tal absorption P = small fc/A turns out to 
be sensitive to the helicity of light: It is substantially 
reduced for one and facilitated for another. Moreover 
changing the excitation energies hio we can investigate 
the dependence shown in Fig. [1] If the out-of-plane 
pseudospin polarization is chosen to be the same in both 
valleys, then the total absorption does not depend on the 
radiation helicity but the two valleys turn out to be differ- 
ently occupied by the photoexcited carriers which is in- 
teresting effect on its own|2^. In the in-plane phase with 
= 7r/2 the total absorption does not depend on light 
polarization, and in the non-interacting limit it equals to 

2 

the universal value as expected [lOj- 

Conclusions. We have demonstrated that the pseu- 
dospin being until now rather uncontrollable and al- 
most unmeasurable quantity can be "unfrozen" by the 
exchange electron-electron interactions ([T]) and play an 
essential role in optical and transport properties of 
graphene. We hasten to say that the Hartree-Fock ap- 
proximation employed here has generically a tendency 
to overestimate ordering such as the pseudospin out-of- 
plane polarization. We believe, however, that the pseu- 
dospin eigenstates Ix^fc) derived above are much more 
robust because their special pseudospin-momentum en- 
tangled structure stems from the free Hamiltonian Hq, 
and the electron-electron interactions do only modify it 



making our predictions reliable at the qualitative level. 
From this point of view the pseudospin can be seen 
as an additional degree of freedom similar to the true 
spin but unaffected by the magnetic field directly. Hav- 
ing this similarity in mind one can think about pseu- 
dospin ferromagnetism ll| , pseudospin accumulation at 
the sample's edge by means of the zero-field Hall current, 
pseudospin selectivity in the optical absorption (|8]), and, 
probably, pseudospin filtering and switching. In a more 
distant future one can imagine some useful effects based 
on the pseudospin polarization like an all-electrical coun- 
terpart for GMR which is obviously very promising for 
application. This Letter should be seen as a first step in 
this direction. 

This work was supported by DFG via GRK 1570. 



[1] K. S. Novoselov et al, Nature 438, 197 (2005). 
[2] A. K. Geim and K. S. Novoselov, Nat. Mat. 6, 183 (2007). 
[3] A. H. Castro Neto et al, Rev. Mod. Phys. 81, 109 (2009). 
[4] C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 

(2005), R. Winkler and U. Ziilicke, Phys. Lett. A 374, 

4003 (2010). 

[5] M. Mecklenburg and B. C. Regan, Phys. Rev. Lett. 106, 

116803 (2011). 
[6] C. Jang et al, Phys. Rev. Lett. 101, 146805 (2008). 
[7] J. Martin et al, Nat. Phys. 4, 144 (2008), E. H. Hwang, 

B. Y.-K. Hu, and S. Das Sarma, Phys. Rev. Lett. 99, 

226801 (2007), S. V. Kusminskiy et al, ibid. 100, 106805 

(2008). 

[8] N. M. R. Peres, F. Guinea, and A. H. Castro Neto, Phys. 

Rev. B 72, 174406 (2005). 
[9] Y. Barlas et al, Phys. Rev. Lett. 98, 236601 (2007). 
[10] J. Nilsson et al, Phys. Rev. B 73, 214418 (2006). 
[11] H. Min et al, Phys. Rev. B 77, 041407 (2008), F. Zhang 
et al, ibid. 81, 041402 (2010), J. Jung, F. Zhang, and 
A. H. MacDonald, ibid. 83, 115408 (2011), J. Jung, et 
al, Phys. Rev. Lett. 106, 156801 (2011). 
[12] N. Nagaosa et al. Rev. Mod. Phys. 82, 1539 (2010). 
[13] N. A. Sinitsyn, J. Phys. Cond. Matt. 20, 023201 (2008). 
[14] N. A. Sinitsyn et al, Phys. Rev. Lett. 97, 106804 (2006). 
[15] N. A. Sinitsyn et al, Phys. Rev. B 75, 045315 (2007). 
[16] W.-K. Tse et al, arxiv: 110 1.2042. 

[17] M. Orlita and M. Potemski, Semicond. Sci. and Tech. 25, 
063001 (2010). 

[18] F. Bonaccorso et al, Nat. Photon. 4, 611 (2010). 

[19] R. R. Nair et al. Science 320, 1308 (2008), A. B. Kuz- 
menko et al, Phys. Rev. Lett. 100, 117401 (2008), K. F. 
Mak et al, ibid. 101, 196405 (2008), E. G. Mishchenko, 
ibid. 103, 246802 (2009), N. M. R. Peres, R. M. Ribeiro, 
and A. H. Castro Neto, ibid. 105, 055501 (2010), F. T. 
Vasko, Phys. Rev. B 82, 245422 (2010), K. F. Mak, 
J. Shan, and T. F. Heinz, arxiv:1012.2922. 

[20] L. O. Juri and P. I. Tamborenea, Phys. Rev. B 77, 233310 
(2008). 

[21] D. V. Khveshchenko, Phys. Rev. Lett. 87, 246802 (2001). 
[22] M. Trushin and J. Schliemann, Europhys. Lett. 83, 17001 

(2008), Phys. Rev. Lett. 99, 216602 (2007) 
[23] D. Xiao et al, Phys. Rev. Lett. 99, 236809 (2007). 



